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Abstract. Temperature and angular momentum dependence of the quadrupole deformation is stud- 
ied in the middle of the sd-shell for 28 Si and 27 Si isotopes using the spherical shell model approach. 
The shell model calculations have been performed using the standard USD interaction and the canon- 
ical partition function constructed from the calculated eigen-solutions. It is shown that the extracted 
average quadrupole moments show a transitional behavior as a function of temperature and the in- 
fered transitional temperature is shown to vary with angular-momentum. The quadrupole deforma- 
tion of the individual eigen-states is also analyzed. 

The study of phase transition in quantum many-body systems is one of frontier research 
topics in physics. Phase transitions are observed both in macroscopic as well as in small 
quantum many-body systems. In macroscopic systems, for instance a metallic supercon- 
ductor, the linear dimension of the system is quite large and the transition from one phase 
to the other occurs at one point [1]. For small systems, the fluctuations play a very cen- 
tral role and the existence of the discontinuity in the heat capacity critically depends on 
the number of constituent particles in the system. This has been demonstrated in small 
metallic grains that discontinuity is observed with large number of electrons in the grain. 
However, as the number of electrons in the grain approaches around 100, the discontinuity 
or the peak structure in the heat capacity disappears [2-5]. 

Phase transitions have also been studied extensively in atomic nuclei using the Hartree- 
Fock-Bogoliubov (HFB) method. The HFB theory predicts phase transition as a function 
of rotational frequency (angular momentum) and temperature (excitation energy). The 
shape transition as a function of rotational frequency has been well studied in most of 
the regions of the nuclear chart. In particular, in the rare-earth region the examples of 
the shape transition are documented in the text books [6,7]. Most of the nuclei in this 
region are prolate deformed with the rotational axis perpendicular to the symmetry axis 
at low spins and it is known that this shape changes to oblate non-collective at higher 
angular momenta in many nuclei in this region. For instance, in the case 160 Yb, the shape 
is prolate for spins up to 40 Ti and above this spin value the shape becomes oblate non- 
collective. The phase transition has also been investigated using the finite temperature 
HFB approach [8,9] and Landau theory [10]. The main conclusion from these studies is 
that nuclei which are deformed at low temperatures exhibit a shape transition to spherical 
shape as the temperature of the system is raised. The critical temperature at which this 
shape transition occurs has a maximum value for isotopes between magic numbers and in 
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rare-earth nuclei this maximum temperature is about 1.85 MeV [10]. 

It is known that the phase transition obtained in the above studies critically depends 
on the inclusion of the quantal and statistical fluctuations [8,11]. In the absence of the 
fluctuations, the phase transition is sharp first order and when some aspects of the fluctua- 
tions are considered, the phase transition is smoothed out. These studies have been mostly 
performed for heavier nuclei using the grand canonical ensemble. The situation is very dif- 
ferent for lighter mass nuclei as it is possible to perform the exact shell model calculations 
and the canonical partition function can be constructed. The canonical partition function 
will incorporate fluctuations accurately. The purpose of the present work is to investigate 
the shape transition in the sd-shell using the canonical partition function constructed from 
the spherical shell model eigen-states. In a recent work, we have used this approach to 
study the pairing correlations [12]. 

The spherical shell model calculations have been performed for 28 Si and 27 Si using the 
recently developed shell model program [12,13]. The shell model Hamiltonian, generally, 
contains single-particle and two-body parts and in the second quantized notation is written 
as 

H = h sp + V 2 , (1) 

where, 

rs 

and 



% = j y~] < rs\v a \tu > c\c\c u c t 
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= V , v / ( 2r + 1 ) = <rs]Va]tu>T 



rstuT 



^4 (rs) x A r (tu^j , (3) 

where e rs are the single-particle energies of the spherical shell model states, which are 
diagonal except in the radial quantum numbers and < rs \v\tu > are the two-body in- 
teraction matrix elements and in the present work are chosen to be those of "USD". The 
two-particle coupled operator in Eq. (3) is given by A\,(rs) = (cjc^)r and Aym t — 
(-l) r ~ Mr A r -M r - The labels r, s, ... in the above equations denote the quantum numbers 
of angular-momentum and isospin of the single-particle states. 'T" quantum number labels 
both angular-momentum and isospin of the two-particle coupled state. The above notation 
is same as that used in [14]. 

In the present work, the statistical averages have been calculated using the canonical en- 
semble approach since the exact solutions have well defined particle number. The average 
value of a physical quantity "F" in canonical ensemble is given by [15,16] 

{{F))=Y J F l e- E ^ T /Z, (4) 
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where, 

=Ei\i), 

Fi = (i\F\i). (5) 

The statistical averages have been evaluated from the lowest one-thousand eigen-states, ob- 
tained from the diagonalization of the shell model Hamiltonian, Eq. (1), for each angular- 
momentum. For a few cases, we have also calculated the averages with 1500 eigen-states 
and results were similar to those with one-thousand eigen-states. 

The average quadrupole moment of the shell model states has been calculated from the 
expectation value of the operator 

Q CpQp &nQn (6) 

where the quadrupole operator is given by 

Q p(n) = v ^r75^(r||r 2 F 2 ||s)(at x a s ) 2 (7) 

rs 

The effective charges e p — 1.5 and e„ = 0.5 have been used and the harmonic oscillator 
length parameter has been calculated with hio — 45A~ 1 / 3 — 25A~ 2 / 3 [17]. 

The shell model calculations have been performed in the middle of the sd-shell for 28 Si 
and 27 Si with the USD interaction [19]. The results for 28 Si and 27 Si are presented in Fig. 
1 with the left side depicting the average quadrupole moment for 28 Si and the right side 
the results of 27 Si. For 28 Si at low temperature, the average quadrupole moment has a 
constant value of about 20 efm 2 up to temperature, T=1.2 MeV and above this temperature 
the quadrupole moment drops. For higher temperatures, the quadrupole moment is noted 
to approach zero and, therefore, indicating that shape transition from deformed to spherical 
shape has occurred. It is noted that the drop in the quadrupole moment is spread from T=l 
to 4 MeV and it is not possible to determine the exact phase transitional point. In the grand 
canonical mean-field calculations, the phase transition occurs at one point. In the present 
canonical study, the phase transition is smeared out due to the fluctuations present. The 
transitional point in quadrupole moment can be approximately inferred as the middle of 
two points - one at which the drop begins and the other point where the constant behavior 
is again observed. These two temperature points for 1=2 in Fig. 1 are 1 MeV and 4 MeV 
and, therefore, the transitional point is inferred as T=2.5 MeV. It is to be noted that positive 
quadrupole moment in Fig. 1 corresponds to the oblate shape for the ground-state band 
of 28 Si using the standard relationship between the laboratory and intrinsic quadrupole 
moments [20]. 

For lighter mass sd-shell nuclei, the question of phase transition has been quite con- 
troversial. The mean-field HFB calculations predict the phase transition, evident from 
the peak structure obtained in the heat capacity [21] at temperature, T=2.1 MeV and the 
vanishing of the quadrupole moment. The quadrupole moment has a finite value at low 
temperatures and then depicts a transition at T=2. 1 MeV. The canonical shell model calcu- 
lations, however, depicted a different behavior. The canonical heat capacity does show a 
peak structure as in the grand-canonical mean-field analysis, but the calculated quadrupole 
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moments didn't vanish. The analysis in these investigations were performed mostly for 
24 Mg and all the angular-momentum states were used to construct the canonical parti- 
tion function [22,18]. In the present study of 28 ' 27 Si, the number of eigen-states for each 
angular-momentum are quite large and the partition function for each angular-momentum 
can be independently constructed. 

It is evident from Fig. 1 that the quadrupole-moments show a transitional behavior for 
each calculated angular-momentum ensemble. The interesting observation from Fig. 1 is 
that the drop in the quadrupole-moment is shifted to lower temperatures with increasing 
angular-momentum. For 1=4, the drop in quadrupole-moment is noted to start at T=0.8 
MeV and the transitional point is inferred to be T=2.3 MeV. For 1=6, the drop in the 
quadrupole moment occurs at a slightly lower temperature as compared to 1=4 and for 1=8 
the drop starts at a very low temperature. We have also analyzed the temperature behavior 
of angular-momentum, 1=1,3, 5 and 7 ensembles and for these ensembles the quadrupole 
moments depict a very irregular dependence on temperature. This would explain the rea- 
son that the quadrupole moments in the earlier study didn't depict a transitional behavior 
as these quantities were averaged over all angular-momentum ensembles. Further, results 
of the present work are in agreement with the mean-field results which predict the transi- 
tional behavior for both quadrupole-moment as well as the heat capacity [18]. Although, 
these mean-field calculations were performed for 24 Mg, but it is expected that these results 
should be similar for 28 Si. 

The results for 27 Si, presented on the right side in Fig. 1, for the angular-momentum 

1 = 5/2, 9/2, 13/2 and 17/2 ensembles, show that the drop in the quadrupole moment has 
a similar behavior as that of the even-even system, 28 Si. The transitional point in the 
quadrupole moment is inferred to be approximately T=2.25 MeV for 1=5/2 . For 1=9/2, 
this transitional point is slightly lowered and occurs at T=2 MeV. For 1=13/2 and 17/2, the 
transition also occurs at about T=2 MeV. 

In Fig. 1, quadrupole deformations have been plotted as a function of temperature for 
a fixed value of angular-momentum. It is also quite instructive to show the angular- 
momentum dependence of the quadrupole moment for a fixed value of temperature. In Fig. 
2, angular-momentum dependence of the quadrupole deformation for temperatures of 0, 1, 

2 and 3 MeV is plotted. It is observed from the figure that as the temperature is raised, the 
quadrupole deformation of higher angular momentum states tends to drop. A clear trend of 
quadrupole moments dropping with increasing temperature for larger angular-momentum 
ensembles is evident from Fig. 2. 

The advantage in the shell model study is that the analysis can be performed for each 
individual state [23]. In Fig. 3, the quadrupole moments are shown for each individual 
eigen-state and have been plotted for the first one-thousand states for both the studied 
nuclei. The quadrupole moment for the low-lying states are randomly distributed from 
large positive to negative deformations. However, with increasing eigen-state number, it is 
evident that the quadrupole moments converge to the spherical shape. In particular, for 1=2 
and 5/2, most of the states at higher excitation energy have zero quadrupole moment. The 
fluctuations from zero quadrupole moment increases with increasing spin and is due to the 
reason that the dimensionality of the basis states becomes progressively smaller with spin. 

In conclusion, temperature and angular momentum dependence of the quadrupole de- 
formation has been studied in the middle of the sd-shell. It is quite evident from the present 
study that collective states for the studied nuclei, 28 Si and 27 Si, depict a transitional behav- 
ior from large oblate deformation to spherical shape. This transition is evident from both 
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canonical ensemble study and the quadrupole moments of the individual states. 

Finally, we would like to mention that the results of the present work are questionable 
due to restricted sd-shell configuration space for T > 2.5 MeV. For higher temperatures, it 
is expected that the fp-shell will be populated. However, the inclusion of fp shell configura- 
tion space is almost impossible in the present context as one needs to evaluate at least 1000 
to 1500 eigen-states for each angular-momentum to have a proper statistical description. 
The shell model Monte Carlo (SMMC) approach [24] is a possible solution to include the 
fp-shell and we intend to look into this problem in the near future. 
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Figure 3. Quadrupole moments are shown for the first one-thousand individual 
eigen-states of the two studied isotopes of 28 Si and 27 Si. 
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